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Abstract: It has been known for some time that a 3D incompressible Euler flow 
that has initially a barely smooth velocity field nonetheless has Lagrangian fluid 
particle trajectories that are analytic in time for at least a finite time (Ph. Serfati 
C. R. Acad. Sci. Serie I 320, 175-180 (1995); A. Shnirelman arXiv:1205.5837 
(2012)). Here an elementary derivation is given, based on Cauchy's form of the 
Euler equations in Lagrangian coordinates. This form implies simple recurrence 
relations among the time- Taylor coefficients of the Lagrangian map, used here 
to derive bounds for the C 1 ' 7 Holder norms of the coefficients and infer temporal 
analyticity of Lagrangian trajectories when the initial velocity is C 1 ' 7 . 
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The issue of finite time blow-up of solutions to the 3D Euler equation for in- 
compressible fluid is still an open question. The available results indicate that (a 
hypothetical) occurrence of singularity is intimately related to the loss of spatial 
regularity (in the Eulerian coordinates) of the solution. It is known, however, 
that an initial regularity of the flow, that is quite modest — marginally better 
than C l , which classical solutions to the Euler equations ought to possess - 
will preserve this regularity for at least a finite time t c (see, e.g., [I])- It may be 
surprising to find that, in such a rough velocity field, the time dependence of the 
position of any fluid particle, i.e. the Lagrangian time dependence, is actually 
analytic up to t c . By contrast, in Eulerian coordinates, lack of spatial smooth- 
ness also translates into lack of temporal smoothness, because of the sweeping 
of fine-scale structure by large-scale velocity fields. For instance, the turbulent 
solutions a la Kolmogorov-Onsager have a spatial regularity of the Holder type 
with an exponent close to 1/3 in both the spatial and the temporal Eulerian 
domain. 
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The unexpected possibility that Lagrangian trajectories could be analytic 
in time was first considered by Serfati [15| (part 2), [T61IT7] (see also [9]) who 
suggested to attack this problem as an abstract ordinary differential equation 
(ODE) in a complex Banach space with an analytic right hand side. Prior to 
this, Chemin j4] proved that the Lagrangian trajectories are C°° smooth. More 
recently, Shnirelman |18) . using the geometric interpretation of the Euler equa- 
tion, obtained another analytic abstract ODE formulation; he then constructed a 
time-analytic Lagrangian solution to the ODE applying Picard's theorem. The 
fact, that the Lagrangian structure can be nicer than the Eulerian one, was 
known to Ebin and Marsden [7] who observed (on p. f04) that in Lagrangian 
coordinates the Euler equations can be written in such a way that "no derivative 
loss occurs" . An early embodyment of this is actually found already in Cauchy's 
little- known Lagrangian formulation of the Euler equation [3J . 

This formulation will be our starting point. It has the advantage of allowing 
to set up not only an elementary proof of the analyticity, but also a procedure 
for calculation of the various temporal Taylor coefficients. In Section^ we derive 
simple recursion relations for the time- Taylor coefficients of the Lagrangian dis- 
placement. Then, in Section|21 we obtain bounds for the C 1 ' 7 of these coefficients. 
The key to obtaining these bound is that the Lagrangian gradient of n-th Tay- 
lor coefficient can be recursively written in terms of nonlocal operators applied 
to products of two and three gradients of lower-order such Taylor coefficients. 
These nonlocal operators are essentially inverse Laplacians composed with two 
suitable space derivatives and thus all the Taylor coefficients stay within the 
same Holder space and simple estimates are obtained for their Holder norms. 
Concluding remarks are presented in Section 2] 

2. Basic equations ... 

Our starting point is not the usual Eulerian formulation of the Euler equa- 
tions for incompressible fluid, but a little-known Lagrangian formulation due to 
Cauchy [3j, specifically his eq. (15), together with the condition of unit Jacobian 
of the Lagrangian map that expresses incompressibility. We use q and x for the 
Lagrangian and Eulerian positions of fluid particles and £ = x — q for the dis- 
placement. In modern notation the dynamical equations of Cauchy [3J take the 
form 
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^V L i fc x V h x k =w , 
fc=l 

det(V L a;) = 1, 

where u>q = V L x Vq is the initial vorticity, V L denotes the gradient in the La- 
grangian variables and V L a; the Jacobian matrix with entries V^Xj. The Cauchy 
equations can be obtained as the (Lagrangian) curl of the Weber equations [TH] 
(see also [5"ll2"irH]V These equations, obviously, express the time invariance of the 
l.h.s. and are therefore sometimes referred to as "Cauchy invariants of an ideal 
fluid" . For simplicity, in what follows we assume spatial 2-7r-periodicity (flow in 
the 3D torus T 3 = [0, 27r] 3 ), but this can be relaxed. 
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Following [5], we study the formal Taylor expansion of the particle displace- 
ment £(q, t) in a power series in time, 

oo 
s=l 



We shall introduce a power scries in t with real non-negative coefficients bound- 
ing certain norms of $} s \q), and prove, under the appropriate assumptions, 
that this series converges for a certain real t — t c . We shall conclude that, for 
complex t, series ([T]) has a finite radius of convergence, not smaller than t c . 

In terms of the displacement, Cauchy's equations become: 
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V L xi + ^V L &xV L a-w , (2) 

k=l 

V L -£+ J2 ((V^)V^-(V^)V^) +det(V L £) = 0. (3) 

l<i<j<3 



Following [8] , we now derive recurrence relations for the time- Taylor coefficients 
of the displacement at t = 0. Substituting JT]) into © and ([3]), we find 

a V L x = u 6t - E «V L e[ n) x V L e^" } , (4) 

fc=l,2,3 
0<n<s 

v l.^)= ^ ((v^v^ 8 -^ - (v t L ei n) )v, L ^ s - n) ) 

l<i<j<3 
0<n<s 

- E ^(V^)^^^" 5 , (5) 

where e^-fc is the unit antisymmetric tensor. It is immediately seen from the two 
equations for s = I that £W = t>o- 

To handle the equations (J4j)— (J5j) for the Lagrangian curl and divergence of the 
Taylor coefficients, we apply the Hodge decomposition, 

£ (s) (<?) = V L x a (s) +V L b {s \ 

where the potentials satisfy the gauge conditions: V L • a/ s ' = and vanishing 
of averages over the periodicity box. This yields a pair of Poisson equations for 
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the vector and scalar potentials and b^: 
sV 2 a (s) = -oj S1 + nyL (( yL x a(n) )k + 

fe=l,2,3 
0<n<s 

x V L ((V L x a^ s ~ n) ) k + Vj^ (s ~ n) ) , (6) 
V 2 6 (s) = (V - ((V L x al">)»+ V^" 1 ) V^(V L x a' 8 -"'),! V^'"" 1 ) 

l<i<j<3 ^ 
0<n<s 

- V s L ((V L x a^) t + V^& (n) )v^((V L x a^-^)j+ Vjb^-^ 
E toiv^^xaW^V^W) (7) 

31,32,33 i=l 
ni+n 2 +n 3 — s 

(here V 2 = (V L ) 2 denotes the Laplacian in the Lagrangian variables and 
(V L x aW)i the fc-th component of the vector V L x a^™'). 

3. ... and their consequence 

In this section we establish the finite-time analyticity in time of the Lagrangian 
map for solutions to the three-dimensional Euler equation. The assumption is 
that the initial velocity field Vq is in C 1,7 . In other words, the initial vorticity 
u>o = V L x Vq is in C ' 7 for some < 7 < 1, i.e., satisfies the Holder condition 
\<*>o{Q + Sq) — u>o(q)\ < C|<5q| 7 . Demonstration of temporal analyticity amounts 
to proving the following statement: 

Theorem. Consider a space-periodic three-dimensional flow of incompressible 
fluid governed by the Euler equation. Suppose the initial vorticity wo(<7) is in 
C 0,7 (T 3 ). Then, at sufficiently small times, the displacement of fluid particles, 
£,(q,t), is given by power series ([I}, whose coefficients can be recursively calcu- 
lated by solving equations (j?]) and ([5]). The radius of convergence of ([T]), T , is 
bounded from below by a strictly positive quantity, which is inversely propor- 
tional to |u> |o, 7 - (More precisely, r > Q c |u;o|o i/2, where Q c is given by (TT2|) 
below.) 

Proof. From recurrence relations ([5]) and ([7]) we derive bounds for suitable 
Holder norms of the potentials. Using our assumption that u)q is in C 0,7 , it is 
easy to establish by induction that each potential aS"' and &W is C 2,7 -regular. For 
this, we recall that C m ' 7 (T 3 ) is an algebra (i.e. the sum and the product of two 
functions from C m ' 7 (T 3 ) belong to this space; in particular, \ fg\o,j < |/|o, 7 |<7|o, 7 ) 
and apply standard inequalities for Holder norms of functions whose integrals 
over the periodicity box vanish: 



\4>U„ < 0|V 2 0| O , 7 V</> G C 2 - 7 (T 3 ), (8) 
|V^V^(<?)|o, 7 < 0|V 2 0| O , 7 and V0 G C 2 < 7 (T 3 ) (9) 
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(see, e.g., [TTl[T2l[T5] or [TO])- Here and are some positive constants depend- 
ing on 7 > and bounded by C/7, V^Vj" denotes the second partial derivative 
d 2 /dqidqj, and |-| s , 7 designates the Holder norm in C S ' 7 (T 3 ) for a scalar function 
4> and |<K<z)|o.7 = max*; \4>k{q)\o.-y for a vector-valued function (p. 
Now, we consider the generating functions 



thus reintroducing time into the problem. Note that, by construction, for t > 0, 
A(t) and B{t) are monotonically increasing functions of time. When no confu- 
sion is possible, we shall write just A and B for A(t) and B(t). In order to show 
that the series A and B converge, we proceed as follows. We derive from rela- 
tions (jSJ) and ([7]) two inequalities for the quantities A and B. The first of these 
inequalities involves time derivatives A and B, but it can be integrated in time. 
The inequality resulting from this integration and the second inequality are of 
the form of upper bounds for A and B. These bounds are polynomials of 2A + B 
that involve, apart from the contribution from the initial vorticity u>o(<z), only 
quadratic and cubic terms. Since and &w vanish at t = 0, the inequalities 
yield an upper bound for 2A + B at least at small times, when 2A + B = O(t) 
and higher-order terms are insignificant. 
Relation © and inequalities imply 



00 



oo 



A(^]T|vV s >|o, 7 i s , 



5(^]T|V 2 &( S >| , 7 ^, 



oc 




s=l fe=l,2,3 
0<n<s 





s>l, 
0<u<s 




\u \ 0n + 66> 2 (2i + B)(2A + B). 



Integration of this inequality in time yields 



A < |wo|o, 7 * + 30 2 (2A + B) 2 . 



Relations ([7]) and inequalities © imply 



B < 60 2 (2A + Bf + 60 3 {2A + B) 3 . 
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From the last two inequalities we deduce 

p(C)=60 3 C 3 + 120 2 C 2 -C + Q>O, (10) 

where we have denoted £ = 2A + B and Q = 2|u>o|o,7 t. The discriminant of the 
polynomial p(£), 

is positive at small times, when Q is sufficiently small, whereby p(Q has three 
real roots The polynomial has then two local extrema at points of different 
signs, hence it has roots of both signs. Since by Viete's theorem the product of 
the roots is negative, two roots are positive and one is negative. Inequality (1101) 
implies 

2A + B<( 2 (Q), (11) 

where C2 is the intermediate root (i.e., the smaller of the two positive roots). 

We determine now the largest time t c , for which bound (jlip holds. Differen- 
tiating the equation for the roots of polynomial (|10p in Q, we find 

dCi (dp 

dQ \d( 

Consequently, on increasing Q, the root £2 monotonically increases and (3 mono- 
tonically decreases till the two roots collide (i.e. (2 = C3) &t a critical value 
Q = Q c . Then t c = Q c /(2|u;o|o,7)- For Q = Q c the discriminant A vanishes, 
which implies 

_2_ 32\ 2 4 2 _2_ 32 

30 + Yj + 276) 2 + 81<9 3 ~ 30 ~ T' ( ' 

For small 7 and hence large ~ 1/7, we have Q c = 6> _2 /48 + 0(6> -3 ) ~ 7 2 . 

So far, inequality pT|) for t < t c = Q c / (2\uj q \ ^) is just an a priori bound. 
However, one can easily check that its derivation can be repeated, without any 
changes, for partial sums of the series A and B having an arbitrary but equal 
number of terms. The bound (|lip being uniform in the number of terms left in 
the partial sums, convergence of the series follows from standard arguments, as 
well as the validity of the bound (TTT|) for the infinite series 2 A + B. This remark 
concludes the proof. 

Clearly, if the bound (fTTj) holds true up to the time t = t c , then the series 
defining A and B are convergent in the complex time disk \t\ < t c . Using the 
inequalities for the coefficients of the series for 2A + B, it is easy to show that, in 
this disk, the series ((TJ) is not just a formal - but a genuine - solution to (EJ)-©- 
By virtue of the proven pQ uniqueness of the solution in the said Holder class, 
any Lagrangian solution takes the form of the series ([1}. This guarantees that, 
within this disk, the solution £ and its gradient V L £ are analytic functions in 
time. 

Finally, we observe that the solution constructed to this point is based on a 
Taylor expansion about t — 0. However, we can restart from any time t > 0, at 
which vorticity has the C ' 7 regularity for some 7 > (which can be smaller than 
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the Holder exponent 7 for the initial data). Thus, temporal analyticity persists 
on the interval [t, t + 1*) for a certain positive i*. This process can be continued 
until the vorticity ceases to have any C ' 7 regularity at some time i», which may 
be finite or infinite. Hence trajectories of fluid particles remain real analytic 
functions of the time as long as the vorticity preserves its Holder continuity. 
This completes our proof of analyticity in time of fluid particle trajectories as 
long as a the velocity field is C 1,7 -continuous for some 7 such that < 7 < 1. 



4. Concluding remarks 

Examination of the structure of the proof of time analyticity given in Section [3] 
reveals that the essential ingredients are (i) the fact that |/g|o.7 < l/lo/ylflio/y 
(algebra property) and (ii) the boundedness in the Holder norm of the operator 
V^V^V -2 , where V~ 2 is the inverse Lagrangian Laplacian acting on functions 
of zero spatial mean. There are obviously many other function spaces for which 
similar properties hold. This will be discussed elsewhere. 

Of course, none of the results derived here from Cauchy's Lagrangian for- 
mulation of the 3D incompressible Euler equation shed light on the issue of 
whether or not there is finite-time blow up. It is interesting to observe that, in 
Cauchy's formulation, all the nonlinear terms vanish if the solutions depend only 
on one spatial coordinate, because all the gradients appearing in the equations 
are then parallel. Our bounds on the various nonlinear terms are rough as they 
do not take into account the cancellation occurring in the vector products, for 
instance, the possible strong depletion of nonlincarity which can be associated 
to the development of almost one-dimensional solutions. At present such issues 
can be investigated only through very precise numerical simulations. By per- 
forming the time-stepping in the Lagrangian domain one can then fully exploit 
the analyticity result. 

Here we have established results for the Lagrangian structure of strong so- 
lutions, as long as they exist. Recently there has been significant progress on 
weak solutions where the velocity has much less spatial smoothness, e.g. Holder 
continuity C ' 7 with 7 close to the Kolmogorov-Onsager value 1/3 [6]. It would 
be interesting to investigate the Lagrangian structure of such weak solutions. 
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